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ABSTRACT 


Terrestrial  spectroscopy  in  the  3-5  pm  wavelength  region  offers  unique  applications,  such  as  detection  of 
carbon  and  sulfur  compounds,  and  subpixel  temperature  determination  of  ground  fires  and  volcanic  lava  flows.  Two 
challenges  in  this  region  are  the  similar  levels  of  solar  reflected  and  terrestrial  emitted  fluxes,  and  the  low  values  of 
both  these  fluxes  making  measurements  of  high  quality  more  difficult  than  in  the  short-wave  or  long-wave  IR 
regions.  We  collected  airborne  hyperspectral  images  with  a  Fourier  transform  spectrometer  over  seawater  in  a 
tropical  environment  combined  with  radiosonde  data  and  MODTRAN  atmospheric  modeling,  to  distinguish  ground 
radiances,  atmospheric  contributions,  and  spurious  instrumental  effects.  Within  the  range  of  our  radiosonde  data  (3 
km  altitude),  MODTRAN  very  closely  reproduces  the  sensor-detected  mid-IR  fluxes.  Clear  separadon  between  solar 
reflected  and  terrestrial  emitted  components  is  achieved,  and  different  atmospheric  composiuonal  spectral  effects  are 
analyzed.  We  applied  the  results  from  this  well-truthed  experiment,  to  spectra  collected  of  the  same  locations  at 
different  times  without  supporting  radiosonde  data,  to  assess  surface  radiance  accuracy  and  atmospheric  variability. 

keywords:  infrared,  spectroscopy,  Fourier,  atmospheric,  MODTRAN 

1  ■  INTRODUCTION 


The  3-5  pm  wavelength  region  is  a  spectral  range  in  the  infi:ared  not  blocked  by  atmospheric  absorptions 
(primarily  from  water  vapor  and  carbon  dioxide).  This  wavelength  region  has  been  used  for  determination  of 
volcanic  lava  flow  and  ground  fire  temperatures,  and  detection  of  carbon,  sulfur,  phosphorous,  and  hydrocarbon 
compounds  15.  A  unique  challenge  for  utilization  of  the  3-5  pm  band  for  remote  sensing  is  the  similar  flux  levels  of 
solar  reflected  and  terrestrial  emitted  radiance.  Unlike  the  visible,  near  infrared  and  far  infirared  where  only  one  or  the 
other  of  these  components  contributes  most  of  the  flux,  in  the  3-5  pm  band  both  components  similarly  contribute  to 
the  total  energy  flux  received  at  a  remote  sensor .  Also,  the  total  radiance  is  low  here,  especially  near  3  pm,  making 
difficult  both  the  detection  and  separation  of  these  radiant  sources  (Figures  1  and  2).  Meanwhile,  as  for  any  remotely 
sensed  measurements  of  the  Earth,  the  atmospheric  transmittance  and  radiance  must  be  accounted  for  to  determine  the 
true  surface  radiant  flux  for  accurate  surface  temperatures  and  compositional  information.  We  have  combined  spectra 
from  the  SMIFTS  Fourier  transform  spectrometer,  temperature  data  collected  in-situ  and  from  low-altitude  radiosonde 
measurements,  and  the  MODTRAN2  radiative  transfer  code.  With  these,  we  have  analyzed  data  in  the  3-5  pm 
wavelength  region  and  applied  it  to  surface  temperature  determination  using  two  methods:  the  first  method  directly 
calculates  an  average  surface  brightness  temperature  from  SMIFTS  sensor-detected  radiance  corrected  using 
MODTRAN2,  and  the  second  method  iteratively  uses  MODTRAN2  with  successive  surface  brightness  temperature 
input  values  to  model  an  at-sensor  radiance  that  best  fits  the  SMIFTS  sensor-detected  spectra.  We  describe  the 
methodology,  and  results  here. 

2.  SMTFTS  INSTRUMENT  AND  DATA  COLLECTION 

We  used  a  spatially  modified  imaging  Fourier  transform  spectrometer  (SMIFTS)®.  Between  3-5  pm,  it  has 
a  signal-to-noise  ratio  of  100-200,  spectral  resolution  of  70  cm'^  and  an  instantaneous  field  of  view  of  0.6  mrad. 
With  SMIFTS,  we  have  collected  mid-wave  infrared  data  over  forest,  desert,  and  tropical  environments. 

Our  data  consists  of  a  series  of  spectral  image  cubes  taken  via  helicopter  approximately  2  km  over  the 
Natural  Energy  Laboratory/Hawaii  Authority  (NELHA)  on  the  island  of  Hawaii.  These  data  were  collected 
throughout  the  daylight  hours  of  May  14,  1994.  NELHA  is  used  for  aquaculture  and  ocean  thermal  energy 


conversion  research,  and  includes  the  maintenance  of  open-air  thermally-controlled  seawater  ponds.  We  used  these 
ponds  as  calibration  for  the  remotely-sensed  helicopter-mounted  SMIFTS  collections,  since  water  is  spectrally  "flat" 
in  the  3-5  pm  region^^.  Average  8-12  pm  pond  brightness  temperatures  were  recorded  by  ground  personnel,  and 
near-ground  air  temperatures  and  relative  humidities  were  also  recorded.  Furthermore,  atmospheric  radiosonde  data 
from  0-3  km  altitude  of  temperature,  pressure,  and  relative  humidity  were  collected  for  one  of  the  pond  overflights. 
Finally,  Global  Positioning  System  (GPS)  location  data  for  the  helicopter’s  latitude  and  longitude  were  available 
throughout  the  day.  With  this  well  instrumented  data  set  (which  may  be  called  "in-situ  assisted  remote  sensing"),  we 
could  focus  on  determining  the  transmittance  and  radiance  of  the  intervening  atmosphere,  to  calculate  true  surface 
flux  from  the  SMIFTS-detected  radiance  for  surface  temperature  comparisons. 

3.  ATMOSPHERIC  AFFECTS  BETWEEN  3-5  urn 

As  described  earlier,  two  main  challenges  for  spectroscopy  between  3-5  pm  are  the  complexity  of  separating 
out  the  solar  reflected  and  terrestrially  emitted  components,  and  the  low  total  signal  from  both  of  these  contributions, 
making  acquisition  of  high  quality  data  and  scientifically  useful  interpretation  difficult.  A  third  concern,  necessary  to 
address  in  the  accurate  determination  of  ground  radiances  and  spectral  features  from  airborne  and  Earth-orbiting 
sensors  in  particular,  is  that  of  correcting  for  the  affects  of  the  atmosphere  between  the  sensor  platform  and  the 
ground  target.  The  solution  to  the  problem  consists  of  characterization,  correction,  and  removal  of  the  affects  of  the 
intervening  atmosphere  on  both  emitted  and  reflected  flux  detected  at  the  sensor. 

Atmospheric  affects  have  only  begun  to  be  quantitatively  addressed  in  the  3-5  pm  spectral  region.  In  the 
visible  to  near-infrared,  several  authors  have  had  good  success.  Blake  and  Singer^  studied  low-level  atmospheric 
correction  using  water  vapor  density  as  a  variable,  and  studies  employing  the  Jet  Propulsion  Laboratory’s  Airborne 
visible/infrared  Imaging  Spectrometer  (AVIRIS)  uses  an  extensive  atmospheric  modeling  process'^’ In  the 
thermal  wavelength  regions,  however,  there  has  been  less  comprehensive  application  of  atmospheric  corrections^, 
and  authors  have  had  to  make  many  simplifying  assumptions  which  were  mathematically  convenient  but  did  not 
yield  accurate  results. 

The  atmosphere  both  attenuates  surface  radiant  flux  and  adds  to  the  total  radiant  flux  measured  by  a  remote 
sensor.  In  the  3-5  pm  spectral  window,  although  the  atmosphere  is  mostly  transparent,  there  are  localized  major 
absorptions  due  to  water  vapor  and  carbon  dioxide,  as  well  as  broad  minor  absorptions  due  to  water  vapor  and  a 
nitrogen  continuum  as  well  as  carbon  monoxide  and  methane,  attenuating  surface  radiant  energy  reaching  a  sensor 
(Figure  3).  Meanwhile,  aerosols  scatter  solar  reflected  and  terrestrial  thermal  energy  as  well  as  emitting  to  add 
radiant  energy.  In  the  3-5  pm  region  for  a  moderate-resolution  spectrometer  at  least  1  km  in  altitude,  not  only  are 
the  at-sensor  solar  reflected  and  terrestrial  emitted  fluxes  similar,  so  are  the  sum  of  those  two  fluxes  and  the  total 
atmospheric  contribution  (Figure  4).  Especially  in  this  spectral  region  of  low  absolute  reflected  and  emitted  solar 
terrestrial  and  relatively  equal  atmospheric  radiances,  it  is  necessary  to  account  for  both  the  atmospheric  extinction 
and  radiance. 

Not  part  of  the  present  focus  but  interesting  to  note,  is  that  the  process  of  removing  atmospheric  affects  to 
determine  correct  ground  radiances  also  yields  useful  information  about  the  intervening  atmosphere  itself,  such  as 
details  about  its  composition  including  water  vapor  and  carbon  dioxide,  and  state  variables  such  as  temperature  and 
pressure.  In  this  way,  a  "sounding"  is  taken  for  the  atmosphere,  as  an  indirect  result  of  accounting  for  its  affects  on 
the  ground  target  spectra. 


4.  MODTRAN  ATMOSPHERIC  MODELING 

Our  primary  tool  to  model  these  atmospheric  affects  was  MODTRAN2  ^ . .  Radiosonde  soundings  were  taken 
once,  in  the  afternoon  of  the  14th,  launched  from  about  7  km  from  NELHA.  The  radiosonde  temperature  and  relative 
humidity  as  a  function  of  altitude  (which  was  calculated  from  pressure)  up  to  3  km  were  input  to  MODTRAN  for 
atmospheric  layer  thicknesses  of  0.2  km,  and  the  MODTRAN  tropical  atmosphere  model  was  used  for  higher 
altitudes.  We  used  the  MODTRAN  multiple  scattering  option:  although  single  scattering  is  sufficient  for  the  most 
part  in  the  mid-infrared,  multiple  scattering  is  important  at  high  solar  zenith  angles  especially  in  the  lower 
atmosphere^,  and  its  inclusion  does  not  appreciably  increase  computer  run  time.  The  MODTRAN  maritime 


boundary  layer  (0-2  km),  spring-summer  (2-30  km),  and  background  "stratospheric”  (>30  km)  aerosol  models  were 
used  without  modification:  there  are  sulfate  (S04^‘)  aerosol  compounds  formed  from  SO2  transported  around  the 
southern  tip  of  Hawaii  and  up  the  western  coast  on  a  160  km  journey  from  degassing  sources  on  Kilauea's  summit 
and  East  Rift  Zone^^,  but  these  compounds  have  only  broad  or  weak  narrow  absorption  features  between  3-5  pm 
overwhelmed  by  the  water  vapor  and  carbon  dioxide  features  Even  if  non-aerosol  SO2  gas  was  dominant,  our 
modeling  shows  that  above-typical  values^  produce  no  measurable  affects  between  3-5  pm  at  70  cm'^  spectral 
resolution.  There  was  an  occasional  thin  fog  haze  blowing  but  no  rain  or  clouds  on  the  1 4th. 

Since  the  radiosonde  balloon  was  launched  some  distance  from  NELHA,  and  because  the  afternoon 
radiosonde  data  was  to  be  applied  at  other  times  of  the  day,  near-ground  air  temperatures  and  relative  humidities  were 
taken  at  NELHA  concurrent  with  each  helicopter  overflight.  The  differences  in  temperature  and  relative  humidity 
between  these  ground  measurements  and  the  surface-layer  radiosonde  values  were  added  to  several  lower  layers  of  the 
radiosonde  inputs  to  MODTRAN,  to  allow  for  the  radiosonde  distance  and  time  differences. 

The  seawater  pond  brightness  temperatures  were  taken  at  NELHA  with  a  Cyclops  hand-held  radiometer, 
concurrent  with  each  helicopter  overflight.  For  the  afternoon  (2:30  P.M.)  data  collection,  which  was  concurrent  also 
with  the  radiosonde  launching,  the  at-surface  brightness  temperature  value  of  one  pond  was  input  to  MODTRAN. 
These  parameters  applied  were  applied  to  MODTRAN,  convolved  to  the  70  cm'^  SMIFTS  spectral  resolution 
(MODTRAN2  highest  spectral  resolution  is  2  cm*^),  and  then  compared  to  afternoon  SMIFTS  radiance  taken  of  the 
same  pond  (Figure  5).  For  the  most  part  the  comparison  of  the  model  run  to  the  measured  radiance  was  very  good, 
convincing  us  that  the  adjustments  made  to  the  lower  atmospheric  layers  of  the  radiosonde  inputs  to  MODTRAN 
were  valid  enough  (at  SMIFTS  spectral  resolution)  to  apply  to  NELHA  pond  overflights  at  the  other  collection 
times  for  surface  temperature  determination. 

5.  SURFACE  TEMPERATURE  DETERMINATION 

The  above-described  MODTRAN  final  inputs,  based  on  afternoon  radiosonde  data  and  known  afternoon 
NELHA  near-ground  surface  temperature  and  relative  humidity,  were  applied  to  helicopter  overflights  with  SMIFTS 
earlier  on  May  14th  to  assess  the  accuracy  of  surface  radiance  measurements  using  our  atmospheric  model 
corrections,  primarily  for  surface  temperature  determination.  The  first  step  in  this  process  was  to  decide  how  best  to 
use  the  MODTRAN-modeled  radiance.  There  are  some  purely  atmospheric  components  to  the  modeled  spectra  which 
contain  no  information  relating  to  surface  radiance,  such  as  the  carbon  dioxide  feature  near  4.3  jim  and  the 
encroachment  towards  5  .um  of  a  water  absorption  band  near  6.7  \im.  There  are  also  some  residual  instrumental 
calibration  artifacts  in  the  measured  spectra,  especially  near  the  increasingly  low  signal-to-noise  towards  3  jiim.  We 
cannot  adjust  the  model  to  account  for  the  instrumental  affects,  but  we  can  account  for  the  atmospheric  affects  by 
excluding  areas  of  the  MODTRAN-modeled  radiance  where  water  vapor  and  carbon  dioxide  absorptions  are  strong. 
For  this,  we  inspected  the  highest  MODTRAN  spectral  resolution  (2  cm'^)  transmittance  (Figure  3)  and  excluded  use 
of  channels  between  4.15-4.5  |im  and  longwards  of  5.25  ixm. 

Using,  then,  regions  of  the  spectra  where  transmittance  is  not  dominated  by  atmospheric  absorptions 
(Figure  6),  we  used  SMIFTS  spectra  taken  of  two  NELHA  seawater  ponds  in  the  early  daylight  hours  of  the  14th 
(6:30  A.M.)  and  applied  MODTRAN  to  determine  their  slightly  different  pond  temperatures  at  that  time.  We  chose 
two  methods  to  do  this,  described  and  compared  below. 

The  first  method  directly  calculates  an  average  surface  brightness  temperature  from  SMIFTS  sensor-detected 
radiance  corrected  using  MODTRAN2,  using  the  following  formula  for  at-sensor  radiance: 

“  (^bb  *  ^w  *  "^a)  +  Lj^  +  Lg  +  Lg 

is  the  total  SMIFTS-measured  at-sensor  radiance, 
is  the  blackbody  emittance  at  the  pond  surface  temperature(s) 
according  to  the  Planck  function, 

is  the  emissivity  of  the  surface,  assumed  0.975  for  seawater^’ 
is  the  MODTRAN-modeled  atmospheric  transmittance  between  the  surface 


where: 


and  the  sensor, 

is  the  MODTRAN-modeled  integrated  atmospheric  path  radiance  at-sensor, 

Lct  is  the  MODTRAN-modeled  ground  reflected  radiance  at-sensor,  and 
L5  is  the  MODTRAN-modeled  path  scattered  radiance  at-sensor. 

The  above  formula  was  solved  for  L55.  This  provided  calculated  radiance  values  for  37  SMIFTS  channels 
between  3.4  -  5.5  jim.  We  then  calculated  the  brightness  temperatures  from  each  channel  radiance  value,  excluding 
those  channels  where  water  vapor  and  carbon  dioxide  absorptions  are  strong.  This  provided  us  with  26  values  of 
surface-equivalent  temperature.  From  these  we  took  an  average.  We  did  this  for  each  of  the  two  seawater  ponds. 

The  second  method  iteratively  uses  MODTRAN2  with  successive  surface  brightness  temperature  input 
values  to  model  an  at-sensor  radiance  that  best  fits  the  SMIFTS  sensor-detected  spectra.  In  other  words,  we  did  not 
use  the  above  formula,  but  completed  the  full  inverse  of  what  we  initially  did  to  match  MODTRAN  to  the  afternoon 
SMIFTS  measurements  from  the  known  afternoon  pond  temperature.  We  took  the  morning  SMIFTS  spectra,  and 
ran  MODTRAN  several  times  at  different  surface  temperatures,  to  find  which  MODTRAN  case  best  matched  the 
morning  SMIFTS  spectra.  Two  best-fit  MODTRAN  input  surface  temperatures  were  found,  one  for  each  pond, 
using  the  method  of  mean  absolute  deviation 

These  methods  of  surface  temperature  calculation  were  also  applied  to  a  helicopter  overflight  of  the  same 
two  ponds  during  late  morning  (1 1:04  A.M.).  Results  from  the  two  methods  for  the  two  ponds  at  the  two  times  are 
shown  in  the  table  below  and  in  Figures  7  and  8.  Also  shown  are  the  "actual’'  pond  temperatures  for  each  time  (in 
Kelvin),  as  measured  by  ground  personnel  with  the  hand-held  radiometer,  as  well  as  the  root-mean-square  deviation 
associated  with  each  temperature. 


Time 

In-situ 

Radiometer 

Temperature 

Method  1 

Method  2 

6:15  A.M. 

296.3  ±  2.4 

296.9  ±  2.7 

297.5  ±  2.8 

6:15  A.M. 

299.7  ±  2.4 

301.4  +  5.1 

301.5  ±  2.8 

11:04  A.M. 

297.3  ±  2.2 

296.9  ±  6.9 

294.2  ±  2.0 

11:04  A.M. 

301.3  ±2.2 

302.3  ±  7.7 

297.7  ±  2.0 

Of  primary  interest  is  that  the  temperatures  determined  using  both  methods  are  consistent  with  the 
radiometer-measured  values,  given  the  rms  deviation  associated  with  each  value. 

Of  secondary  interest  is  to  note  that  method  1  gives  the  overall  closest  determined  temperatures  to  the 
radiometer-measured  values,  even  though  method  2  gives  the  overall  least  root-mean-square  deviation  in  each  of  those 
temperatures.  These  points  follow  from  the  nature  of  each  calculation.  Method  1  is  calculated  from  the  actual 
SMIFTS-detected  radiance  -  which  is  more  representative  on  average  of  true  surface  temperature  even  though  it 
contains  channel-specific  instrumental  and  atmospheric  "noise,"  Method  2  is  calculated  from  MODTRAN-modeled 
radiance  -  which  although  it  contains  no  "noise"  and  so  produces  a  smoother  radiance  without  spurious  variation,  is 
still  modeled  data  and  so  is  likely  to  be  less  accurate  on  average  than  measured  data. 

It  is  useful  to  compare  the  temperatures  and  rms  deviations  determined  by  each  method  to  present 
capabilities  of  the  NOAA  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  aboard  the  NOAA  9  satellite. 
That  instrument  is  a  standard  yardstick  for  terrestrial  surface  temperature  determination  in  the  infrared.  Kalluri  and 
Dubayah^  describe  three  methods  of  surface  temperature  calculation  from  AVHRR  data,  giving  rms  values  between 
AVHRR  and  ground  measurements  of  2.01,  2.20,  and  3.60  K,  which  are  consistent  with  the  generally  accepted  ±3  K 
absolute  accuracy  range  for  AVHRR.  Higher  accuracies  have  been  obtained  with  AVHRR^  but  these  have  been 


in  tandem  with  data  from  other  sensors  and/or  numeric  algorithms  which  are  not  yet  commonly  accepted.  Therefore, 
compared  to  standard  AVHRR  surface  temperature  rms  accuracies,  method  1  above  ranges  from  at  least  as  good  to 
somewhat  worse  than  AVHRR,  while  method  2  fares  at  least  as  good.  It  would  be  difficult  to  compare  the  absolute 
temperature  determination  accuracies  of  the  above  two  methods  to  AVHRR  without  a  more  complete  analysis  of  the 
accuracy  of  our  hand-held  radiometer  and  the  instrumentation  employed  for  the  AVHRR  ground  truthing. 

6.  SUMMARY  AND  CONCLUSIONS 

Recent  advances  in  technological  and  computational  capabilities,  combined  with  the  use  of  the  most  current 
version  of  radiative  atmospheric  modeling,  is  allowing  useful  spectroscopy  in  the  3-5  Lim  infrared  region.  This 
region  was  previously  avoided  due  to  low  absolute  levels  of  total  radiation  flux  and  the  difficulty  of  separating  the 
solar  reflected  and  terrestrial  emitted  fluxes.  Our  initial  measurements  using  SMIFTS  combined  with  MODTRAN2 
atmospheric  modeling  have  resulted  in  remotely  sensed  sea  water  temperature  determination  consistent  with  ground- 
based  measurements,  and  within  the  same  accuracy  ranges.  Finally,  with  further  improved  technology  and  use  of 
MODTRAN3  soon  to  be  available,  not  only  will  more  accurate  remotely  sensed  surface  radiances  be  calculable  in 
this  spectral  window,  but  very  useful  3-5  jim  atmospheric  absorption  bands  such  as  from  carbon  and  sulfur 
compounds  will  be  utilizable  for  purely  atmospheric  applications. 
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9.  FIGURES 


FIGURE  1.  Representative  terrestrial  at-surface  radiant  fluxes  (without  specific  atmospheric  absorption 
features),  showing  the  relative  energy  levels  of  reflection  and  emission. 
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FIGURE  2.  Representative  terrestrial  at-surface  radiant  fluxes  (without  specific  atmospheric  absorption 
features)  in  between  3-5  jam,  showing  the  low  flux  level  of  both  components,  as  well 
as  their  relative  equivalence. 
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FIGURE  3.  M0DTRAN2  2  cm**  spectral  resolution  representation  of  atmospheric  transmittance 
showing  molecular  absorptions  primarily  due  to  water  vapor,  gaseous  carbon  and 
nitrogen  compounds. 
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FIGURE  4.  Comparison  of  at-sensor  solar  reflection  (100%  albedo  surface),  terrestrial  (blackbody) 
emission,  and  atmospheric  emission  (all  at-sensor).  Although  our  surface  was  seawater 
(assumed  2,5%  albedo),  the  above  illustrates  the  general  equivalence  of  the  three 
components  contributing  to  the  radiance  in  this  wavelength  region. 
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FIGURE  5.  MODTRAN-modeled  and  SMIFTS -detec ted  spectra  for  2:30  P.M.  May  14,  1994, 

approximately  2  km  over  NELHA  seawater  pond  of  known  surface  temperature. 
The  MODTRAN  run  included  the  use  of  concurrent  radiosondes,  and  adjustments 
to  several  lower  atmospheric  layers  of  radiosonde  inputs  to  account  for  differences 
in  temperature  and  relative  humidity  between  NELHA  and  the  radiosonde  launch. 
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FIGURE  6.  Figure  5  repeated  to  show  regions  of  the  modeled  and  measured  spectra  where  the 
atmospheric  transmittance  is  not  dominated  by  atmospheric  absorptions,  used 
for  surface  temperature  determination. 
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6:15  A.M.  (top)  and  11:04  A.M.  (bottom)  pond  temperatures  via  ground 
measurement  (hand-held  radiometer),  effective  blackbody  calculation  (method  1), 
and  reverse-MODTRAN  modeling  (method  2),  and  the  associated  root-mean-square 
errors.  Note  the  y  axis  scale  offset. 
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